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We investigate the superfluid state of repulsively interacting three-component (color) fermionic 
atoms in optical lattices. When the anisotropy of the three repulsive interactions is strong, atoms 
of two of the three colors form Cooper pairs and atoms of the third color remain a Fermi liquid. 
An effective attractive interaction is induced by density fluctuations of the third-color atoms. This 
superfluid state is stable against changes in filling close to half filling. We determine the phase 
diagrams in terms of temperature, filling, and the anisotropy of the repulsive interactions. 
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The cold atoms in optical lattices offer a new way 
of studying fascinating aspects of quantum many-body 
physics [l-yl . High controllability of the atoms allows us 
to create novel systems, which extends our current knowl- 
edge in, e.g., condensed matter physics. The research 
on Bose-Fermi mixtures trapped in optical lattices ex- 
tends the paradigm of the Mott transition Q ■ Fermionic 
atoms with odd numbers of internal degrees of freedom 
are other examples. Recently, interesting features of de- 
generate three-component fermionic gases were actively 
investigated 0-0] • For three-component (color) fermionic 
atoms in optical lattices, it was shown theoretically that 
the paired Mott insulator (PMI) can appear in spite of 
half filling, where the average atom number per site is the 
non-integer 3/2 [1, When two of the three repulsions 
between atoms with different colors are stronger than 
the other, pairs of weakly repulsing atoms are formed to 
avoid the two stronger repulsions. Therefore, in the PMI 
the effective particle number at a site becomes an integer 
Q. This feature implies that, in the Fermi liquid (FL) 
states in the vicinity of the PMI transition point at half 
filling and close to half filling, superfluid pair fluctuations 
are strongly enhanced at low temperatures, although the 
bare interactions are repulsive. We expect a character- 
istic superfluid (SF) state to appear there. In this SF, 
atoms of two of the three colors form Cooper pairs and 
atoms of the third-color remain a FL. This SF state is 
similar to the color superfluid of the attractively inter- 
acting three-component systems [Iol - [l2j |. Knowledge of 
this novel SF provides an insight into the Cooper pair- 
ing mechanism of other repulsively interacting systems 
in, for instance, condensed matter physics. 

In this Letter, we investigate the SF transition in re- 
pulsively interacting three-component fermionic atoms in 
optical lattices. We first discuss the effective interactions 
using Feynman diagrammatic approaches. By combining 
the dynamical mean-field theory (DMFT) with a pertur- 
bative approach [3] , we then investigate the SF proper- 
ties in detail. We show that the SF phase appears close to 



the PMI phase, when the anisotropy of the interactions is 
strong. We reveal that the effective attractive interaction 
that forms Cooper pairs of two-color atoms is mediated 
by density fluctuations of the third-color atoms. We also 
perform unperturbative calculations using the self-energy 
functional approach [14| . The two sets of results are com- 
plementary. 

The low-energy properties of the present system are 
well described by the following Hubbard Hamiltonian: 
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where the subscript is the summation over the 

nearest- neighbor sites, and a ia (di a ) and hi a are creation 
(annihilation) and number operators of a fermion with 
color a at the ith site. Here, t denotes the hopping in- 
tegral. We choose the chemical potential [i a to satisfy 
the balanced population of each color atoms; n{= (n a )). 
For simplicity, we set U12 = U and U23 = t/31 = U', and 
Ml = A*2 = A* an d M3 = We introduce an anisotropy 
ratio R = U/U'. We neglect the confinement potential 
for the first approximation. We also neglect the symme- 
try breaking phases except for the superfluid state. 

Let us first consider the effective interaction (U12) be- 
tween color- 1 and 2 atoms by adopting a random phase 
approximation analysis as shown in Fig. [1] The contri- 
bution of all the bubble diagrams is summarized in the 
following expression, 
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FIG. 1: (Color Online) Diagrammatic representation of the 
effective interaction between color-1 and 2 atoms. 



where x = [ivi, k), vi is the bosonic Matsubara frequency, 
and k is the wave vector. The bubble diagram corre- 
sponding to density fluctuations is given by x«( x ) = 
- J2 y 9a(y)g a {y + x), where y = (iu n , q) with w„ and q 
being the fermionic Matsubara frequency and the wave 
vector, respectively, and the bare Green's function of 
color-a atoms is denoted by g a (y) = l/(iw n + A*q — £q)- 
The other effective interactions are evaluated in the same 
way. Since U\ 2 = U and U\$ — U23 = U' , we set 



X2 = X an d X3 = x' ■ The effective interactions 



U(x)(= C/i 2 (x)) and C/'(x)(= t/ 23 (x) = C/ 3 i(x)) are given 
by 



L>(x) 



«7-C/'V(x) 



[1 - C/x(x)][l + C/x(x) ~ 2l/' 2 x(x)x'(x)] 

tr 

l + C/x(x)-2C/' 2 x(x)x'(x)' 



(4) 



For {/ -C U' < t (R -C 1), U(x) ~ C7 - £/'Y(x). The 
second term can be regarded as the effective attraction. 
When the effective attraction overcomes the bare inter- 
action U, t/(x) becomes attractive. By contrast, for 
[/={/'< t (R = 1) we obtain the effective repulsive 
interactions t/(x) = f7'(x) ~ U. 

This simple discussion suggests that when the 
anisotropy is strong (R <C 1) the effective attractive in- 
teraction between color-1 and 2 atoms is induced by den- 
sity fluctuations of the color-3 atoms. Furthermore, for 
R <C 1 the FL-PMI transition occurs in the strongly cor- 
related region as shown in Ref. 8]. Therefore, we can 
expect the s-wave SF caused by the attractive £/(x) to 
appear in the FL phase close to the PMI transition point 
at low temperatures. We perform DMFT calculations to 
confirm this expectation. When the system approaches 
the PMI transition point, local correlation effects en- 
hance the quantum fluctuations. In particular, SF pair 
fluctuations between color-1 and 2 atoms are expected 
to be strongly enhanced. To deal precisely with the lo- 
cal correlation effects we apply the infinite-dimensional 
system, which is a nontrivial limit for discussing both 
the Mott transition and the s-wave superfluidity 1J, 15| • 
For noninteracting atoms, we use a semicircular density 
of states, p (x) = yftt 2 - x 2 /{2nt 2 ). 

To investigate the SF, we carry out the second-order 
perturbative expansion of the self-energy within the 
DMFT framework. This theory is namely the iterated 
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FIG. 2: (Color Online) Diagrams for the (a) normal and (b) 
anomalous self-energies of color-1 and 2 atoms, respectively, 
and (c) the normal self-energies of color-3 atoms. 



perturbation theory (IPT) [la |. By using the Nambu 
spinor representation ip q = (c q i, cLq 2 i c q3)i the Green's 
function matrix is given by 



Gi(y) F(y) 
G(y) = I F(y) -Gi(-y) 

G 3 (y) 



(5) 



In the infinite-dimensional system, the self-energy matrix 
is purely local and independent of q, yielding S(y) = 
H(iu) n ), where 

(Ei(iw„) S(iu n ) \ 

S(icj n ) -Ei(-io; n ) (6) 
S 3 ( iWn ) / 

According to the DMFT [13], S(iw„) can be ob- 
tained from the effective impurity model, and the im- 
purity Green's function matrix Gi mp (iui n ) must be self- 
consistcntly equal to the local Green's function matrix of 
the original system, Gi oc (iui n )[= I] q G(y)]. 

We calculate the first- and second-order self-energies, 
£W and S^ 2 \ respectively, by collecting the diagrams 
shown in Fig. [2] Here, the arrowed lines denote the 
diagonal (g ca ) and off-diagonal (/ c ) components of the 
cavity Green's function matrix g c , which is the nonin- 
teracting impurity Green's function matrix and satisfies 
the Dyson equation, g c — [Gr + S] _1 . To obtain pre- 
cise results close to the PMI transition point, we modify 
the self-energy matrix as S = I]' 1 ) + AS^ 2 ' according to 
Ref. [IB} . The matrix A is determined so as to reproduce 
the exact 1/u) expansion of Gi mp . The diagonal elements 
are A n = A 22 = [R 2 n 2 (l - n 2 ) + n 3 (l - n 3 ) - i? 2 |$| 2 + 
2 J R(Z) 23 -n 2 n 3 )]/[i? 2 n c2 (l-n c2 )+n c3 (l-n c3 )-i? 2 |$ c | 2 ], 
^33 = [ni(l-?^i)+£>i 2 -?iin 2 +|$| 2 ]/[7i cl (l-n c i)+|$ c | 2 ], 
and the off-diagonal elements are zero. The number of 
atoms n a and the SF order parameter $ = (fijia^) are 
calculated from the diagonal and off-diagonal elements of 
Gi mp (iu} n )e luJn0 , respectively. In the same way, n ca 
and <I> C are calculated from the cavity Green's function. 
The double occupancy D a p = {n a np) is obtained from 
D Q/3 = -Tr[G imp dE/dt/ a/3 ]. 

By using the IPT, we calculate <&, D a p, and the renor- 
malization factor Z a = 1/(1 — dT, a /dui)\ LJ= o. Figure[3]^a) 
and (b) show the U'/t dependence of these quantities at 
R = 0.1 and T/t = 0.03 for different fillings n = 0.5 
(half filling) and n — 0.48, respectively. In Fig. [3£a), 
Z a decreases with increasing in U'/t and vanishes at 
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FIG. 3: (Color Online) (a) and (b) Renormalization factor Z a , 
SF order parameter <E>, and double occupancy D a p calculated 
with the IPT as functions of U'/t at R = 0.1 and T/t = 0.03 
for fillings (a) n = 0.5 and (b) n = 0.48. (c) Eigenvalues (k) 
of K ab at T/t = 0.03 (top) and 0.02 (bottom) for R = 0.1 
and n = 0.48. Negative values denote the instability of the 
homogeneous phases. For comparison, $ is also plotted. 




FIG. 4: (Color Online) (a) Anomalous self-energy ReS(w+), 
ImS(aj+), and imaginary part of the polarization function 
Imx3(kJ+). (b) Single-particle excitation spectra p a (oj). 
These quantities are calculated with the IPT. Parameters are 
U'/t = 3, R = 0.1, n = 0.48, and T/t = 0.01. 



U'/t ~ 4.2. The double occupancy D\i (-D13) increases 
(decreases) as U'/t increases. These results indicate that 
the PMI transition occurs to avoid the stronger U'/t. 
When U'/t approaches the PMI transition point from be- 
low, the SF order parameter becomes finite and reaches 
zero with a tiny jump at the transition point, suggest- 
ing a discontinuous transition. In contrast, $ increases 
smoothly at U'/t ~ 2.5 with increasing U'/t. The tran- 
sition between the FL and SF is thus continuous. For 
n = 0.48 the PMI does not appear as shown in Fig. [3jb) , 
because the effective particle number at a site is non- 
integer [1, [lf|. Note that the density- wave states are 
also unstable, unless at least two of three color atoms are 
at commensurate filling 0, Gjl EH- On the other hand, 
the SF phase appears in 3 < U'/t < 4.2 for n = 0.48. 
We thus expect this characteristic SF to be stable near 
half-filling for the balanced population. 

We further examine the stability of the above SF phase 
against the phase separation (PS), which was pointed out 
in recent theoretical studies on the attractively interact- 
ing three-component systems We consider the 
possibility of the PS between paired color-1 and 2 atoms, 
and unpaired color-3 atoms. To this end, we calculate 
the 2x2 compressibility matrix K a b = dN a /dpb with 
N a = ni + ri2,n3 and p,b = /i, p! . Figure [3Jc) exhibits 
the eigenvalues of K a b at T/t = 0.02 and 0.03 for R = 0.1 



and n = 0.48. We find that for both T/t = 0.02 and 0.03 
one of the eigenvalues shows anomaly at U'/t ~ 4.8 and 
becomes negative in U'/t > 4.8, which indicates the in- 
stability to the PS. In contrast, the stable homogeneous 
SF and FL phases are found in U'/t < 4.8. 

We now investigate the SF in terms of its dynami- 
cal properties. As discussed concerning Eq. ([3]), density 
fluctuations of color-3 atoms play a key role in the ap- 
pearance of the SF. Therefore, we discuss Imxs(w) and 
the dynamics of the anomalous self-energy S(u>). In Fig. 
|4ja) , we show Imx3(w+), ReS'(aj + ), and Im5(cj+) with 
lu + = lu + i0+ for U'/t = 3, R = 0.1, n = 0.48, and 
T/t = 0.01. The spectra Imx3(w) have a peak around 
ujp/t ~ 1.6, indicating that the scattering caused by the 
density fluctuations of color-3 atoms occurs prominently 
around to ~ uj p . The spectra ReS'(aj + ) exhibit resonant 
behavior and change their sign at co/t ~ 1.8, which is 
larger than oj p /t. At u/t ~ 1.8, lmS(ui+) has a peak. 
Since ReS'(u;+) is a measure of the effective interaction 
21], the effective interaction between color-1 and 2 atoms 



is attractive when uj/t < 1.8. The frequency of the peak 
of ReS'(cj + ) is smaller than oj p /t, and ReS'(cj + ) remains 
negative for a large to ft. These features are caused by 
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the bare repulsion U between color-1 and 2 atoms 
These dynamical properties are qualitatively the same as 
those of the strong-coupling superconductor caused by 
the electron-phonon interaction 



21] 



In Fig. |4jb), we show single-particle excitation spectra 
defined as p a (io) = —(l/w)Ixn.Gi mPta (uj+). The s-wave 
superfluid gap is seen around the Fermi energy (lu = 0) 
in p\(uj)[= P2 (w)], while the renormalized peak is seen 
in ps(uj) owing to the FL nature. This s-wave spectral 
gap is consistent with the dynamical properties shown 
in Fig. HJa), which are qualitatively the same as those 
of the strong-coupling phonon-mediated superconductor. 
We find that the incoherent Hubbard band structures ap- 
pear in both spectra px{oS) and p^(uj), which stem from 
strong correlation effects. Actually, the peak position 
of lm.S(u!+) is in this incoherent Hubbard band region 
of Pi(uj). On the basis of the results described so far, 
we can argue that the effective attractive interaction be- 
tween color-1 and 2 atoms is mediated by the density 
fluctuations of color-3 atoms. 

By systematically changing T and U' , we determine 
the finite-temperature phase diagrams for n — 0.5 and 
0.48 at R = 0.1. The comparison of these two phase 
diagrams is very informative, although we neglect the 
density-wave state at n = 0.5. For n = 0.5, the SF phase 
appears adjacent to the PMI phase as shown in Fig. [5][a). 
For n — 0.48, the SF phase extends to a stronger U'/t 
region than that for n = 0.5 as shown in Fig. [5jb), sug- 
gesting that SF pair fluctuations survive in the vicinity 
of the PMI phase. However, the homogeneous phase be- 
comes unstable in U'/t > 4.8. The inset in Fig. Mjo) 
shows the temperature dependence of <& for n — 0.48 
and R = 0.1. We find $ cx y/T -T c with T c being the 
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FIG. 5: (Color Online) Phase diagrams obtained with the 
IPT. U —T phase dia grams for i? = 0.1 and different fillings 
(a) n — 0.5 and (b) 0.48. Inset in (b) shows $ as a function 
of T/t at 7? = 0.1 for n = 0.48. (c) Phase diagrams for 
T/t — 0.03. Main panel and inset are U'-R and U'-n phase 
diagrams for n = 0.48 and R — 0.1, respectively. 
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FIG. 6: (Color Online) (a) SF order parameter <I>, Double 
occupancy D a p, (b) renormalization factor Z a , and (c) grand 
potentials of the SF and the normal states (FL and PMI). 
These quantities are calculated using the self-energy func- 
tional approach for 7? = 0.1, T = 0, and n = 0.5. 



SF transition temperature, which is typical of weak- and 
strong-coupling phonon- mediated superconductors [2l|. 
We also determine the R- and n-dependent phase dia- 
grams shown in Fig. G3c). From Eq. ([3]) we expect 4> 
to decrease and then vanish with increasing i?, because 
the effective attractive interaction cannot overcome the 
bare repulsion U in the weak anisotropy (large R) region. 
We find that the SF appears in R < 0.11 as shown in 
the main panel of Fig. 03c), suggesting that the strong 
anisotropy is thus necessary for the appearance of the 
SF. We remark that as n deviates from 0.5, <E> decreases 
monotonically due to the suppression of density fluctu- 
ations of color-3 atoms. We find that the SF survives 
down to n ~ 0.43 as shown in the inset of Fig. [5Jc). 

Next, we discuss the SF transition using an unper- 
turbative numerical method. The higher-order contri- 
butions, which are not included in our diagrammatic 
analyses, may play an important role. We thus per- 
form the self-energy functional approach calculations 
[l4l | . This method was successfully used for studying the 
attractively interacting three-component Hubbard model 
[TlL H3 |. In this method we introduce a proper refer- 
ence system, which has to include the same interacting 
Hamiltonian as that of the original Hubbard Hamilto- 
nian. Wc employ the two-site Anderson impurity model 
as the reference system. The local correlation effects 
are unperturbatively taken into account via the refer- 
ence self-energy calculated by the exact diagonalization 
method [U, [l| [3 . In Fig. Ufa) and (b), we show $, 
D aj3l and Z a for n = 0.5, R = 0.1, and T = 0. The 
U'/t dependences of these quantities are qualitatively the 
same as those obtained by the IPT. In Fig. EJc), we com- 
pare the grand potentials f2 of the SF and the normal 
state obtained by the calculations without any symme- 
try breaking. Note that the normal state in U'/t < 3.2 is 
the FL, while that in U'/t > 3.2 is the PMI. The grand 
potential of the SF is lower than that of the FL. We thus 
confirm that the SF is stable. The qualitative agreement 
between the IPT results and the present unperturbative 
calculations means that our IPT treatment captures the 
essentials of the correlation effects. 



Finally, we discuss the relation between our theoreti- 
cal results and the experiments. The double occupancy 
shown in Fig. |3ja) and (b), and the single-particle ex- 
citation spectra shown in Fig. Eflb) can be measured by 
state-of-the-art experimental techniques 0, HH. Atoms 
with convenient magnetic Feschbach resonances such as 
6 Li arc good candidates for realizing the SF, because 
the anisotropy of the interactions is controllable. It was 
pointed out that, for attractively interacting 6 Li atoms 
in optical lattices, the strong three-body loss makes the 
homogeneous phases unstable fioj ] . However, as shown in 
the experiment [5j, the loss rate in the repulsive region 
is significantly small as compared with that in the at- 
tractive region. Furthermore, the three-body loss effects 
are suppressed in the repulsively interacting system, be- 
cause the two-body repulsions decrease the triple occu- 
pancy (nxriins). Accordingly, the loss-induced PS can 
be negligible and the phase diagrams shown in Fig. [5] are 
adequate for the repulsively interacting three-component 
6 Li atoms in optical lattices. 

In summary, we have pointed out that the SF appears 
in repulsively interacting three-component fermionic 
atoms in optical lattices. We have shown that the 
Cooper-pairing mechanism is basically similar to that 
of the conventional phonon-mediated superconductivity, 
while the SF phase appears close to the PMI phase. 
When we consider the realistic situations, we have to 
take account of the confinement potential and the effects 
of the small but finite loss rate depending on species of 
atoms. These issues are subjects for a future study. 
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